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We present molecular dynamics study of a generic (coarse-grained) model for 
single-polymer diffusion confined in a corrugated cylinder. For a narrow tube, 
i.e., diameter of the cylinder 5 < 2.3, the axial diffusion coefficient Dn scales as 
Du oc iV~ 3 / 2 , with chain length N, up to N ~ 100 then crosses over to Rouse scal- 
ing for the larger N values. The A r_3//2 scaling is due to the large fluctuation of 
the polymer chain along its fully stretched equilibrium conformation. The stronger 
scaling, namely iV -3 / 2 , is not observed for an atomistically smooth tube and/or for 
a cylinder with larger diameter. 



The dynamics of polymer chains in high geometric confinements has attracted increased 
attention over the last decades l|-|3|]. While early studies striven to achieve a fundamental 
understanding of the often exotic and unpredictable behavior of polymers in confinement, 
recent research has mostly been motivated by applications that range from biology to tri- 
bology 4j-t7j. Despite significant progress within the field, in particular on the structural 
and thermodynamic properties of polymers in confinement, their dynamics have remained 



poorly understood. In particular, the studies are mostly done 



mers 



14( and/or polymers confined between two surfaces 15l4l8l|. In this context, both 



or strongly adsorbed poly- 



the chain length dependent mobility 8|, [9|, as well as the effect of varying surface coverage 
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concentration are studied 



101 ] . Furthermore, the studies of these systems are not only re- 



DNA, proteins and phospholipids |8|, lUl 



stricted to the synthetic polymers |9l, [10|, b ut also to biologically relevant systems, such as 



Despite significant interest in studying the polymer dynamics in confinement, a lesser 
investigated system is the polymers in cylindrical confinement. In this context, one of 
the pioneering works to study the statics and dynamics of polymer confined in a cylinder 
employed Monte-Carlo simulations |l9 |. Studies also include the investigation of static prop- 
erties, such as swollen to globule transition, in a soft tube 20] and the relaxation dynamics 
of polymer chains confined in a cylinder 21] . In a more recent work, it was shown that the 
DNA confined in a nanotube can exhibit large fluctuations at the shorter length scales and 
these fluctuations are ceased at longer length scales 22|. Furthermore, as far as the scaling 
laws of the chain length dependent dynamics is concerned, axial dynamics in a confining tube 
is usually said to follow Rouse dynamics, i.e., diffusion coefficient D\\ oc iV -1 with the degree 
of polymerization N 2^. In this context, the computational studies have usually employed 
ideally flat confining surfaces. However, what has not (yet) been investigated is the extent 
to which the corrugation or surface roughness affects the polymer dynamics in a cylindrical 
confinement. Corrugation, more generally speaking breaking translational invariance, is a 
necessary "ingredient" to exert shear force during the polymer dynamics in high geometric 
confinements. In this context, simulati ng gen eric bead-spring polymers moving on surfaces 



that only had atomic-scale corrugation [12J, ll3[ , it was shown that the model could reproduce 
perhaps the counter-intuitive experimental work of lateral dynamics of polyethylene glycol 
(PEG) on fused silica surface by Zhao and Granick 10], in which the lateral dynamics D 
first increased with the surface coverage concentration V and then suddenly dropped to a 
small value at a threshold concentration V* . The simulations suggested that this change in 
D may be due to a structural transition from single to double layers and that the double 
layers have more geometric flexibility to lock into the substrate's registry, which increases 
energy barriers and thus reduces lateral mobility at large V {12I \\\. 



In a previous Letter [141 ]. one of us have shown that the lateral diffusion D of a strongly 
adsorbed polymer chain onto a corrugated surface gets strongly influenced by the surface 
roughness. The simulations showed D oc N~ 3 / 2 with chain length N, if the polymer is 
incommensurate with the substrate which is consistent with the experiments of polymer 



adsorbed on a solid surface 



9j. Furthermore, Rouse scaling could be found if the polymer 
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is commensurate with the substrate and/or for a smooth surface, as observed for polymer 
dynamics on a soft surface [8|. The weaker scaling is a result of vanishing friction correlation 
for a polymer adsorbed on a commensurate or on a smooth surface. 

Motivated by the success of our previously used model Jjj], which could reproduce 



the generic features of the existing experiments [SHlO], we revisit the system of polymer 
dynamics under cylindrical confinement and study the effect of atomic level corrugation on 
axial diffusion Dm. More specifically, we will investigate the scaling of Dn in the presence of 
surface roughness. 



In this work, we use the well known bead-spring model [241 ] . The individual monomer 
units interact via a truncated repulsive Lennard- Jones (LJ) potential, with a short range 
cut-off r c = 2 1//6 . This r c is chosen to model good solvent condition as well as potential 
is continuous at r c = 2 1 / 6 . Results are presented in terms of the LJ length scale a, the 
LJ energy parameter e and the mass of individual monomer m. Adjacent monomers in 
the polymer chain interact via an additional Finitely Extensible Nonlinear Elastic (FENE) 
potential, 

Wr)4-^ ln|1 - ( ^ for ^ (1) 
J oo elsewhere 

where Rq = 1.5 and k = 30. The parameters of the potential (i.e., LJ + FENE) ensure 
no unphysical bond crossing is allowed and give rise to an effective bond length of around 
0.97. The bead-spring polymer chain is confined in a tube that has atomic scale roughness. 
For this purpose we choose couple of surface architectures. In one case, we choose a surface 
that has the same inherent structure as a carbon nanotube and in the second case we roll 
a triangular lattice into a cylinder. The surface structures and the corresponding cylinders 
are presented in the Fig. [TJ Periodic boundary conditions are applied only along the axis of 
the tube, i.e., z direction. Monomer and tube atoms interact via a repulsive LJ potential 
with a cutoff 2 1 / 6 . The sizes of the polymer chains are varied from N = 20 to N = 500. 
We choose several different tube diameters 5 in both the systems (we will provide details at 
the appropriate captions). For every N, the length of the tube is chosen at-least twice the 
size of the fully stretched polymer conformation. To avoid commensurability effects [14] , the 
nearest neighbor atomic distance is chosen to be 0.8 for carbon nanotube architecture and 
1.209 for triangular lattice. To mimic static free energy barrier we keep the surface atoms 
frozen. Moreover, it is yet important to mention that we have also conducted a few test 
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FIG. 1: (color online) Snapshot of the surface architecture used in the simulations, part (a) repre- 
sents a triangular lattice and part (b) is a replica carbon nanotube. 

runs where the surface atoms were coupled to their respective lattice sites with a harmonic 
spring with a spring constant k = 5000. In this case also the results were found to be the 
same as the case of rigid surface atoms. 

Temperature is imposed through a Langevin thermostat. Since the substrate atoms 
are constrained to their equilibrium positions, the motion of the polymers needs to be 
thermostated so that constant temperature conditions can be mimicked. For this purpose, 
we chose a Langevin thermostat, which only acts on the motion normal to the axis of the 
cylinder, i.e., normal to the polymer axial dynamics. The coupling strength is chosen as 
7 = 0.001, which ensures a diffusive behavior within the choices of 5. This value of 7 
is chosen after calculating diffusion coefficients for different 7's and making sure that the 
diffusion coefficient does not get affected by the specific choice of 7 (data not shown). Unless 
stated otherwise the thermal energy is set to 0.5. It is yet important to mention that, since 
our calculations are performed in thermal equilibrium and thus in linear response, damping 
due to wall friction and damping due to hydrodynamic interactions are linearly additive. 
Therefore, under extreme confinement, as in our case, the surface damping will dominate 
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and will screen out the hydro dynamic damping. Hence, we abstain to include explicit 
solvent in order to emphasize more on the wall friction. However, in some cases, influence of 
hydrodynamic interactions can be visible on the lateral dynamics of polymers near smooth 
surfaces, where hydrodynamic scaling, namely iV _3//4 , can be observed [251] . Furthermore, as 
expected, hydrodynamic interactions are screened for a rough surface |25[. In this study, the 
configurations are equilibrated for a few 10 7 MD time steps (depending on the chain size) 
and then observations are carried out over another 20 x 10 6 MD time steps; this corresponds 
to a time of 10 5 in the LJ units. 
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FIG. 2: (color online) Diffusion coefficient Dm along the tube function of chain length N. 

Part (a) shows the data for the cylinder with carbon nanotube architecture for two different tube 
diameters 5. The thermal energy is elevated to 0.7 for 5 = 2.1. We show the results for the default 
system (rough surface) as well as the result for a smooth surface where the LJ radius between 
monomer and surface atom interaction was increased to 1.5. In the case of smooth surface, we 
choose 5 = 2.6 to have the same accessible are for the monomer as the case of rough surface and 
5 = 2.1. In part (b), we show the results for cylinder with triangular structure. Diameter was 
chosen as 5 = 2.3 and the thermal energy is elevated to 0.7. Note that the data sets in two plots 
are obtained from two different (independent) molecular dynamics softwares. For part (a), we used 



DL-POLY 



261 ] and for the part (b) we used a self-written parallel MD code. 
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The main goal of this paper is to study the dynamics of polymers under cylindrical 
confinement. However, before describing dynamics we also want to ensure that our model 
reproduces correct scaling for the static properties. A quantity that best describe the struc- 
tural property is the radius of gyration R g (data not shown). We find the radius of gyration 
along the axis follow R g u oc N. Furthermore, for 5 < 2.3 the monomer fluctuations normal 
to the tube axis are almost negligible, which ensures an almost one dimensional polymer 
conformation, giving a hint for the single file motion in a tube. 

Now we focus on the polymer dynamics, in particular on the axial diffusion coefficients 
D\\. For this purpose, we start with the calculation of the mean square displacement of the 
individual monomers C\(t), 

1 N 

ci(t) = -j2(m+t , )-R i (t')] 2 ), (2) 

i=l 

and the center-of-mass C^), 

C 2 (t) = (jR{t + t>) - R{t')] 2 ) , (3) 

where Ri is the coordinate of the i th monomer along the tube and R is the center-of-mass 
coordinate of the polymer chain along the tube. The lateral diffusion coefficient is related to 
the C(t), which shows a linear time dependence at larger times, i.e., 2dD = lim^oo dC(t)/dt 
with d as the dimensionality. For monomer mean square displacement, we observe an initial 
ballistic region where Ci(i) oc t 2 when t < 1, then between 1 < t < 10 2 the system exhibits a 
sub-diffusive behavior with C(t) oc t 1 / 2 . Finally, for t > 10 2 , the system goes to the diffusive 
region. This behavior is observed for N = 100. However, for larger N values, diffusive region 
starts at a bit longer time scales. On the other hand, Ci(t) reaches the diffusive region much 
earlier than the monomers and directly crosses over from the ballistic to the diffusive regime. 
For all the chain lengths N, the D» was calculated by measuring the slope of C^t) at large 
times. 

In the Fig. [21 we summarize the most important result for D\\ for both systems. In part 
(a) we present the results for replica carbon nanotube and in part (b) we have the system 
with triangular geometry. The data for 5 = 2.8 in part (a) of Fig. [2] are consistent with 
the scaling law D oc iV -1 predicted by Rouse model and with prior simulation results j^ . 
It is also important to mention that dynamics of polymer between two surfaces usually 
supports the Rouse scaling [27]. However, upon decreasing the inter- wall distance to a value 
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where the polymer conformation is ideally flat, one can expect a stronger scaling, namely 
D oc iV _3//2 , similar to that of the scaling law obtained for the strongly adsorbed polymers 
[jjj]. Though, by no means, we want to make any comparison of cylindrical confinements 
with the adsorbed and/or confined polymer cases. The apparent deviation from the Rouse 
scaling invokes the question of whether a iV -3 / 2 scaling can also be observed in the case of 
a polymer confined in a corrugated cylinder. 

In order to investigate the effect of corrugation and possible observation of stronger 
scaling, the diameter of the confining tube is reduced to a value of 5 = 2.1. As indicated 
earlier, the polymer conformation is almost one dimensional for 5 = 2.1, similar to that of 
single file motion. We also want to point out that due to the "physical" presence of the 
surface atoms, actual accessible space for the monomers normal to the surface is around 
unity in L J units, which is also the size of the monomers. Indeed a scaling law of D oc iV~ 3//2 
is observed for N < 100. However, for N > 100, interesting enough, we see an "unexpected" 
crossover to D oc N^ 1 (see Fig. EJ^a)). Furthermore, in order to check the correctness of 
this cross-over dynamics, we run independent simulations on another surface architecture 
where the surface topography was chosen to be of triangular lattice. In Fig. [2](b), we show 
the results for Du as a function of N. Here we choose 5 = 2.3. Indeed, a cross-over from 
D oc iV~ 3 / 2 to Rouse scaling is again observed. Furthermore, we see that cross-over occurs 
at around N = 100. (Note: two distinctly different softwares are used to obtain these two 
sets of simulation data). While the cross-over scalings are always of great interest, it is yet 
imperative to understand the molecular origin of the behavior, which is directly accessible 
in molecular dynamics simulations. As of the iV -1 scaling at large values of N, an expected 
argument would be that polymers may move as uncorrelated domains and thus lead to Rouse 
dynamics. Since a polymer chain is in a single file conformation, the polymer is thus rubbing 
against the tube where almost every monomer is in contact with the tube wall, so that the 
damping (i.e., the inverse diffusion constant) of the polymers centroid is simply proportional 
to the number of monomers in contact with the cylinder. Therefore, for strong confinements 
the argument naturally leads to Rouse dynamics. 

Now we look into the iV~ 3//2 scaling: first intuition would suggest that the motion in 
this region is correlated and therefore leads to a stronger scaling. Indeed while looking at 
the simulation snapshots, we see a fluctuating motion of the center-of-mass of the polymer 
around its stretched equilibrium configuration along the tube axis (similar to that of a 
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FIG. 3: (color online) Normalized probability distribution P{5z) of the axial hopping distance Sz 
of center-of-mass after every 10 LJ time units for N = 50 and 100 LJ time units for N = 500, 
respectively. 

trajectory of the periodic motion in a moving vehicle as observed from a stationary reference 
frame). On the other hand the motion is rather smooth for N > 100, as expected by the 
observation of Rouse motion in this regime. 

To make a more quantitative illustration we have drawn Fig. |3l where we show the 
probability distribution P(Sz) of the axial hopping distance 5z of center-of-mass after every 
10 LJ time units for N = 50 and 100 LJ time units for N = 500, respectively. It can 
be appreciated that the significant hopping is observed in the centroid motion for N = 
50, which is visible from the long tail in the probability distribution function. It is yet 
important to mention that the fluctuations can be as big as 2 to 3 lattice distances. Whereas 
the fluctuations are rather small and/or minimal for the large iV's given the large time 
interval over which the hopping distance is calculated. These fluctuations enhance the overall 
correlation along the chain and hence leading to N~ 3 / 2 scaling law in the mobility. We have 
also measured the end-to-end relaxation r as a function of chain length. As expected, r oc iV 3 
for Rouse scaling and r oc N 3 - 5 for D oc iV~ 3 / 2 . 

We also want to emphasize that the fluctuations are observed due to the coupling of 
transverse and axial motions in high geometric confinement, which are "only" due to the 
presence of atomic scale roughness. This affects the dissipation of kinetic energy associated 
with the polymers center-of-mass motion along the axis. Therefore, if our scenario is correct, 
it would be important to validate if the said iV~ 3 / 2 scaling would disappear by employing 
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an atomistically flat cylinder and/or under moderate confinement. As of the later case, we 
have already shown that D oc iV _1 for 5 = 2.8 (see Fig. EJ^a)). To mimic the smooth surface 
we have increased the LJ radius between monomer and surface to 1.5. We also choose larger 
S that ensures same normal accessible space for the monomer in the smooth tube as that in 
the case of the rough tube. As shown in Fig. EJ^b), data clearly reveal a iV -1 scaling law. 
Also, as expected, fluctuations are not observed for the smooth surfaces. 

In conclusion, using molecular dynamics simulations of a generic model, we have studied 
the dynamics of polymer chain trapped in a "replica" carbon tube and a corrugated tube 
with triangular geometry. We introduced atomic scale roughness to take account of the 
friction between polymer and cylindrical tube. We observe a cross-over scaling law for the 
axial dynamics Du of the polymer chain as a function of chain length N, namely N~ 3 ^ 2 to 
N" 1 . We argue that the cross-over happens due to the surface corrugation that leads to 
the coupling of tranverse and axial motions, leading to large axial fluctuations of the center- 
of-mass for chain lengths iV < 100. Whereas, the fluctuations are ceased for the longer 
chains and hence the Rouse scaling is observed. The cross-over disappears for the tube 
diameter longer than the monomer size and/or for a smooth surface. The results presented 
in this work might be of particular importance for the study of the translocation of polymers 



through a narrow channel. In this case, the "so called" scaling behavior of the trans 
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time t will strongly depend on the dynamics within the cylindrical confinement 
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